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The dramatic increase in the efficiency of a quan- 
tum computer over a classical computer 1,2 , raises 
a natural question asking, how much of this suc- 
cess could be attributed to its quantum nature 
and how much to its probabilistic content. To 
highlight this issue, we put forward the novel 
idea of a possible chemical computer driven by 
reaction-diffusion (RD) processes 3-5 based on a 
probabilistic but classical approach. Such com- 
puters, obeying non-equilibrium statistical me- 
chanics, can describe superpositions of empty 
and filled states with certain probabilities. With 
these probit states serving as computational ba- 
sis states, such RD computers with operations 
satisfying a necessary semi-group property could 
mimic some well known quantum logic gates and 
carry out teleportation like procedure using en- 
tangled states, believed to be a prerogative of the 
quantum world. Moreover, assuming a nonlinear 
extension the RD computers could be used for 
cloning of arbitrary states, which is a famous for- 
bidden operation in standard quantum computa- 
tion. 

It is important to compare the performance of a quan- 
tum computer (QC) with that of a classical but stochas- 
tic computer, to understand more clearly the reason be- 
hind the success of a QC. Though the randomised algo- 
rithm of Solovay-Strassen has prompted generalisation of 
the universal Touring machine to a probabilistic one 1,2 , 
apparently not much attention has been paid in suggest- 
ing a computer, driven by probabilistic processes, as a 
possible alternative to a QC. We observe however that, 
such classical stochastic events, could be realised in phys- 
ical systems controlled by RD processes 3-5 , which can 
describe a state as the superposition \P >= po\0 > 
+Pi|l >, with pa + pi = 1, where po is the probabil- 
ity with which an empty state |0 > (o) can occur, while 



Pi is the corresponding probability of |1 >, i.e of the 
state occupied by a single particle (•). Such processes 
may be involved in an one dimensional array of N lattice 
sites and along with diffusion to the right or to the left 
(i + ohoIi), can describe reactions like coagulation, 
decoagulation (• + •<->•) , birth and death (o <-> • + •), 
occurring with certain probabilities. The evolution of 
these probabilities is governed by a master equation 
dtPcr = Ho-tPt, where p CT is the probability of the system 
to be in state \cr = (ci, 02 . . . , <jn) >, defined by a con- 
figuration stretched over N lattice sites with er^'s being 
or 1. The dynamical operator i? CT1 - expressed through 
transition rates between the configurations er and r is as- 
sumed to be a linear operator determined by the physical 
processes involved. By introducing \P >= X) CT P CT I <T - > ' 
the evolution equation can be cast as a real Schrodinger 
equation dt\P >= H\P >, with the state evolving in a 
finite time as \P{t) >= U(t - t )\P(t ) >. The linear 
evolution operator U with only positive and real matrix 
elements has an important constraint 
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imposed by the conservation of total probability: 
J2aP<? = 1- We call any matrix satisfying condition 
as the RD matrix. 

We can notice the striking similarity between these re- 
lations and those of the quantum mechanics, where only 
the superposition coefficients are the probability ampli- 
tudes, the evolution equation is the complex Schrodinger 
equation and the evolution operator is a unitary ma- 
trix. This similarity, though formal, should be motiva- 
tion enough to explore the idea of a RD computer for a 
possible alternative to a QC. 

Note that, for an operation to qualify for any se- 
quential computation, a primary requirement is that, it 
must exhibit at least a semi-group property. The quan- 
tum computation is guaranteed to have the semi-group 
property, following from its unitary group transformation 
linked to the reversibility of quantum dynamics. Return- 
ing to the RD processes, we find rather unexpectedly 
that, the defining constraint Q is preserved under evo- 
lution: EiUijit-to) = £ fc (Ei(M*-*i))tffcj(*i-*o) = 
J2kUkj{ti — to) = 1 and since U = I satisfies 
the set of RD matrices defining RD operations forms 
a semi-group. Moreover, since the inverse operation 
(Z7) -1 , if it exists, should also satisfy QJ, some of the 
RD operations might form a group. However, in gen- 
eral the required positivity condition on the matrix ele- 
ments are not preserved under inverse operation and the 
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group structure is obtained formally, only if this positiv- 
ity condition is relaxed. We find that, such a group with 
N x N matrices with real entries satisfying (JJJ would 
involve N(N — 1) generators and act as a transforma- 
tion group on a plane embedded in a iV-dimensional 
Euclidean space, with defining relations J2i $i = 1 
on its coordinates. Such groups are simplest examples 
of algebraic groups. Interestingly, the matrices satisfy- 
ing additional constraints similar to Q on their rows 
: ^2f Uij = 1, i = 1, ■ • • , N constitute a subgroup of di- 
mension (iV-l) 2 , which acts as a stationary group keep- 
ing the point with x\ — x 2 — ■ ■ ■ — xn — on the plane 
fixed under its action. 

Therefore, we may apply in principle the RD processes 
for computation purposes and in analogy with a qubit 
call the elementary information in a RD computer, based 
on probabilistic approaches, as probit. In contrast to the 
qubits belonging to a Hilbert space with the state vec- 
tors describing a unit sphere, the probit states would 
belong to a real vector space, with the tip of the state 
vectors moving on a plane limited between positive coor- 
dinate axes of unit lengths. The action of RD matrices 
keeps this plane invariant by transforming a point to an- 
other point on the same plane. For single probit states 
the RD operators are 2x2 matrices defined by two in- 
dependent parameters t/oo = a, Uoi = b. We find a re- 
markable property that hm„_ >00 ([/'(a, b)) n — U(a oc , a^), 
where = 1 _ b a+b . That is on repeated operation by an 
arbitrary RD matrix, all state vectors end up in an equi- 
librium state \s >= OoojO > +(1 — aoo)|l > ; irrespective 
of the probit state we start with (see Fig 1). 




Figure 1 Single probit states under action of RD operations, confirming 
also the semi-group property of RD matrices. Under repeated operations by 
U(a, b) all vectors from arbitrary initial position j-q evolve into a fixed point 
r g — (aoo, 1 — Qoa) (see text). For symmetric matrices, this fixed point 
r s = (i, ^) becomes independent of even the operating RD matrix. In the 
figure particular values a = 0.4,6 = 0.2, r = (0.7, 0.3) , r g = (0.25,0.75), 



and for the symmetric case a = 0.4, b = 0.6 have been considered. 

For demonstrating the action of such computers we 
start with simplest RD operations, which could act like 
logic gates and transform computational probit states 
>,|1 > to mixed states: |0 >^ U O o\0 > +i7io|l > 
, |1 >^ Un\l > +U O i\0 >, with U being a RD matrix 
satisfying Q). Notice that, if Uoi — t/io = 1, an empty 
state ends up giving birth to a particle, while a filled 
state becomes empty through a death process, i.e. the 
operation acts like a NOT-gate, similar to a quantum 
logic (QL) X-gate. Therefore, an arbitrary single pro- 
bit state \P >= po\0 > +pi\l > after passing through 
such a NOT-gate, due to the linearity of RD operations 
would yield \P >= pi\0 > +po\l >■ i- e the probabil- 
ities of computational states would toggle. Note that, 
other popular QL gates like Z-g&te and the Hadamard- 
gate having negative matrix elements 1 are not allowed 
here and this particular limitation, as we see below, pre- 
vents a RD computer to match a QC in solving some 
important algorithmic problems. Reciprocally however, 
while only the NOT-gate and the identity operation can 
have equivalent QL gates, all other values of Uij would 
represent new types of single probit gates achievable by 
RD computers alone. For example, the matrix with all 
Uij = |, having no inverse, transforms both |0 >, |1 > 
to the same uniform state \u >= |(|0 > +|1 >). We use 
below this uniforming RD gate in designing teleportation 
of probit states. 

Focusing now on two-probit operations we show that, 
C-NOT and Swap gates can also be defined through RD 
computation similar to QC. Defining in analogy with QL 
gates 1 the first probit as the control state \x > and the 
second one as the target state \y >, we may construct a 
C-NOT gate through RD operations, such that, if \x = 
> is an empty state, nothing happens to \y >, but if a 
particle is present in \x = 1 >, the computational state 
\y > is negated through the NOT gate presented above. 
We check easily that the matrix jj CNOT representing this 
operation satisfies the RD condition QJ. Observe also 
that, the control states can be copied by fixing y = 0, 
only if they are in the computational basis and not in 
the arbitrary state \P >. This indicates that the no- 
cloning theorem, a well known no-go theorem in QC, is 
valid also for the RD computers. However we will see 
below that, allowing nonlinear RD processes one could 
bypass this theorem. 

We may construct the swap logic gate P driven by 
a RD computer, again in analogy with the QL gate. 
The action of such a gate through RD operations would 
mean, when the neighbouring two sites of the system are 
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the same, i.e, either both empty or both filled- nothing 
changes. However, when a particle appears only at one 
site, it diffuses to the empty site, either to the right or 
to the left. We check again that this is a permitted RD 
matrix operation satisfying Q. 

More interesting f-gates 1 can be constructed as well in 
a RD computer, where a function x — > f(x) mapping the 
values 0, 1 to themselves, is induced by a logic gate, that 
changes the target state \y >^ \y+f(x) >, depending on 
the nature of the control state \x >. We see immediately 
that similar to the C-NOT gate all such f-gates satisfy the 
RD condition (JTJ. Note that, for all above logic gates like 
NOT, C-NOT, Swap and f-gates the related RD matrices 
have well-defined inverses and therefore each of these RD 
operations forms a group. Moreover such RD operations 
formally coincide with the corresponding QL gates ex- 
hibiting unitarity and therefore the RD computers using 
these gates can serve as possible alternative to QC for 
related computations. 

Another operation where RD computation can mimic 
a QC is the celebrated teleportation problem 1 . Here 
also the probabilistic logic can go parallel to the quan- 
tum logic and as in the quantum case we may consider 
that, Alice and Bob share an entangled two-probit Bell- 
like state |(|00 > +|11 >). For teleporting a RD state 
\P >= po\0 > +pi\l >, Alice passes it as a control state 
through a C-NOT gate together with the first probits of 
the shared entangled state in her possession as the tar- 
get, similar to the procedure adopted in QC. This gives 
the three probit state ±(p (|000 > +|011 >) +Pi(|110 > 
-4- 1 101 >)), the first probits of which are then sent by 
Alice through a uniforming RD gate, elaborated above, 
yielding the final state §{(|00 > +|10 >)\P > -+-([11 > 
+ |01 >)\P >} with \P > being the toggled state. There- 
fore, if now Alice in her measurement of her two probit 
states finds either |00 > or 1 10 > the state \P > itself 
has been teleported, while if the result is either 1 1 1 > 
or 1 01 >, the toggled-state \P > has been sent, in which 
case the original state is recovered by passing it through a 
NOT-gate. The whole process could also be repeated us- 
ing another Bell-like state |(|01 > +|10 >), which is the 
only other entangled state possible in RD computation. 
Thus it appears that, quantum teleportation like scenario 
could be reproduced within the classical RD processes. 

Limitations of RD computers would show up in deal- 
ing certain algorithmic problems, which renders efficient 
solutions in QC. Deutsch problem is to tell whether a 
Boolean function / is constant (i.e. /(0) = /(l) ) or bal- 
anced ( i.e. /(0) 7^ /(l)) by testing the outcome of an 
operation, while its Jozsa generalisation asks the same 



question (with slightly different definitions) for a N-bit 
state 1 . For testing these problems in a RD computer 
we may take the uniform state \u > as the control and 
|0 > as the target and see that, this pure input state after 
passing through the f-gate remains pure (either \u > |0 > 
or \u > |1 >) if the gate is constant, but turns into an en- 
tangled state in the balanced case. Therefore if somehow 
we can test for the entanglement, that would solve our 
problem. However, how to do this efficiently is not clear 
at this stage. Similarly, Grover's search algorithm 1,2 in 
RD computation may also start with the uniform state 
\u > passing through the oracle (initialised in the state 
|0 >). If the oracle is attached to a search function / 
yielding f(X) — 1 for the solution set X and 0, other- 
wise, one naturally gets at the output a superposition 
of orthogonal states J2x \ X > I 1 > + J2x' \ X ' > 1° >■ 
Therefore, though the solution states \X > are distin- 
guishable theoretically from the other states \X' >, ap- 
parently in no way we can detect them to match the effi- 
ciency of a QC. In case of Shor's factorisation algorithm 1 , 
the case is even worse, since from the very beginning it 
uses complex superposition coefficients, which are not 
permitted in RD operations. 

In spite of these deficiencies of the proposed RD com- 
puter, its domain permissible beyond reversible processes 
might give some generalities, as we have seen with the 
uniforming RD gate. Another possibility of extending 
RD processes with important potential applications is to 
induce nonlinear operations. Let us suppose that we in- 
tend to clone an arbitrary state \P >= p \0 > +pi\l > . 
Since we have the indication that this is not possible 
with linear RD processes, we allow the RD operator 
U = I <g) U(P) to depend on the initial state \P >. If 
it acts on the product state ^\P > (|0 > +|1 >), sim- 
ple calculations reveal that, one could obtain the cloned 
state \P > \P > at the output point, by requiring the 
nonlinear dependence as /7(P)oo + U(P)qi = 2p . and 
U(P) W + U(P) ii = 2p v We notice that the RD ma- 
trix condition still holds due to po + p\ — 1. Note 
that, since in RD processes the transition rates are de- 
termined mostly phenomenologically, the assumption of 
nonlinearity should not lead to any serious consequences, 
like violation of causality, faster than light propagation 
etc. as one must expect with similar nonlinear assump- 
tion in quantum computations 1 . 

Therefore similar to quantum computers, our pro- 
posed classical computers based on reaction-diffusion 
processes and using probabilistic probit states can con- 
struct many important logic gates and carry out telepor- 
tation like processes, while through their nonlinear exten- 



3 



sion the cloning of mixed states could be made possible. 
Though the present work is purely theoretical, this novel 
concept of RD computation, being classical in nature, 
might be easier than a quantum computer to realise and 
handle in practice. 
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